Statistical mechanics approach to the probability distribution of money 
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This Chapter reviews statistical models for the probability distribution of money developed in 
the econophysics literature since the late 1990s. In these models, economic transactions are 
modeled as random transfers of money between the agents in payment for goods and services. 
Starting from the initially equal distribution of money, the system spontaneously develops 
a highly unequal distribution of money analogous to the Boltzmann-Gibbs distribution of 
energy in physics. Boundary conditions are crucial for achieving a stationary distribution. When 
debt is permitted, it destabilizes the system, unless some sort of limit is imposed on maximal debt. 

"Money, it's a gas." Pink Floyd, Dark Side of the Moon 



I. INTRODUCTION 

In this Chapter, the probabihty distribution of money 
in a system of economic agents is studied using 
methods of statistical physics. This study originates 
from the interdiscip linary field known as econophysics 
( Stanley et all Il996l) , which applies mathematical meth- 
ods of statistic al physics to so cial, economical, and finan- 
cial problems (|StauffeiJ . [2003) . 

One of the most puzzling social problems is the persis- 
tent wide range of economic inequality among the pop- 
ulation in any society. In statistical physics, it is very 
well known that identical ("equal") molecules in a gas 
spontaneously develop a widely unequal distribution of 
energies as a result of random energy transfers in col- 
lisions between the molecules. Using similar principles, 
this Chapter shows how a very unequal probability distri- 
bution of money among economic agents develops spon- 
taneously as a result of money transfers between them. 

The literat ure on s o cial a nd economic inequality 
is enormous ( Kakwani 19801). Man y pape rs in the 
economic lite r ature (jChampernownel . Il953l : ICibratl . 
I1931I : iKaleckil Il945() use a stochastic process to de- 
scribe dynamics of individual wealth or income and 
to derive their probability distributions. One might 
call this a one-body approach, because wealth and 
income fluctuations are considered independently for 
each economic agent. Inspired by Boltzmann's kinetic 
theory of collisions in gases, we introduce an alternative, 
two-body approach, where agents perform pairwise 
economic transactions and transfer money from one 
agent to another. We sta rt with a simple pairwise- 
transf er model proposed by iDragulescu and Yakovenkg 
(|200Cll ). This model is the most closely related to the 
traditional statistical mechanics, which we briefly review 
first. Then we discuss other money-transfer models 
and further developments. For a more detailed and 
systematic review of the progress in this field, see 
the recent review by Yakovenko and Rosscrl (|2009|), a s 



well as reviews by Chatteriee and Chakrabartil 
iRichmond. Hutzler. Coelho. and Repetowica 



2007 ): 



2006 ) 



[Richmond. Repetowicz. Hutzler. and Coelhol (|2006[ ) and 
a popular article bv iHavea (|2002| ). 

Interestingly, the study of pairwise money transfer and 
the resulting statistical distribution of money have vir- 
tually no counterpart in the moder n economic literature . 
Only the search theory of money ( Kivotaki and Wrightl . 
I1993I ) is somewhat related to it. However, a probabil- 
ity distribution of money among the agents within the 
search-theoretical approach was onl y recent l y obt ained 
numerically by the economist Miguel iMolicd ( 20061 ). 



II. THE BOLTZMANN-GIBBS DISTRIBUTION OF 
ENERGY 

The fundamental law of equilibrium statistical me- 
chanics is the Boltzmann-Gibbs distribution. It states 
that the probability P{e) of finding a physical system or 
subsystem in a state with the energy e is given by the 
exponential function 
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where T is the temperature, and c is a normalizing con- 
stant ([Wannier . 1987i ) . Here we set the Boltzmann con- 
stant ks to unity by choosing the energy units for mea- 
suring the physical temperature T. Then, the expecta- 
tion value of any physical variable x can be obtained as 
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where the sum is taken over all states of the system. 
Temperature is equal to the average energy per particle: 
T ^ (e), up to a numerical coefficient of the or der of 1. 



Eg . (Ill) can be derived in different ways ( Wannien . 
Il987| ). All derivations involve the two main ingredients: 



statistical character of the system and conservation of 
energy e. One of the shortest derivations can be sum- 
marized as follows. Let us divide the system into two 
(generally unequal) parts. Then, the total energy is the 
sum of the parts: e = £1 + £2, whereas the probability 
is the product of probabilities: P{e) = P{£i) P{£2)- The 



only solution of these two equations is the exponential 
function ([1]). 

A more sophisticated derivation, proposed by Boltz- 
mann, uses the concept of entropy. Let us consider N 
particles with the total energy E. Let us divide the en- 
ergy axis into small intervals (bins) of width Ae and 
count the number of particles Nk having the energies 
from Ek to Ek + Ae. The ratio N}./N = Pj, gives the 
probability for a particle to have the energy Sk- Let us 
now calculate the multiplicity W, which is the number 
of permutations of the particles between different energy 
bins such that the occupation numbers of the bins do 
not change. This quantity is given by the combinatorial 
formula in terms of the factorials 



Am to the agent j for some goods or services, the money 
balances of the agents change as follows 



W 



Nl 
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The logarithm of multiplicity is called the entropy S — 
In W. In the limit of large numbers, the entropy per 
particle can be written in the following form using the 
Stirling approximation for the factorials 
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^PfclnPfc. (4) 



Now we would like to find what distribution of particles 
among different energy states has the highest entropy, 
i.e., the highest multiplicity, provided the total energy of 
the system, E = J^k ^k£k, has a fixed value. Solution of 
this problem can be easily obtained using the method of 
Lagrange multipliers (jWannier , 1987 ) , and the answer is 
given by the exponential distribution ([T|). 

The same result can be also derived from the er- 
godic theory, which says that the many-body system oc- 
cupies all possible states of a given total energy with 
equal probabilities . The n it is straightforward to show 
( Lopez- Ruiz et all . l2008l ) that the probability distribu- 
tion of the energy of an individual particle is given by 
Eq. (P). 



III. CONSERVATION OF MONEY 

The derivations outlined in Sec.|TT]are very general and 
only use the statistical character of the system and the 
conservation of energy. So, one may expect that the ex- 
ponential Boltzmann-Gibbs distribution ([T]) would apply 
to other statistical systems with a conserved quantity. 

The economy is a big statistical system with millions 
of participating agents, so it is a promising target for ap- 
plications of statistical mech anics. Is there a conserved 
quantity in the economy? iDragulescu and Yakovenkd 
POOO) argued that such a conserved quantity is money 
m. Indeed, the ordinary economic agents can only re- 
ceive money from and give money to other agents. They 
are not permitted to "manufacture" money, e.g., to print 
dollar bills. Let us consider an economic transaction be- 
tween agents i and j. When the agent i pays money 



m.j — > m'i = ruj — Am, 



rrij — > 171;= rrij 



(5) 



The total amount of money of the two agents before and 
after transaction remains the same 

rrii + rrij = m[+ m'j , (6) 

i.e., there is a local conservation law for money. The rule 
([5]) for the transfer of money is analogous to the trans- 
fer of energy from one molecule to another in molecular 
collisions in a gas, and Eq. ^ is analogous to conser- 
vation of energy in such collisions. Conservative models 
of this kind are also studied in some economic literature 
(iKivotaki and Wrightl flQOl iMolicol 120061 . 

We should emphasize that, in the model of 
IDragulescu and Yakovenko' [200u} [as in the economic 
models of [Kiyotaki and Wright (1993); Molico (2006)], 
the transfer of money from one agent to another repre- 
sents payment for goods and services in a market econ- 
omy. However, the model of IDragulescu and Yakovenkd 
(2000) only keeps track of money flow, but does not keep 
track of what goods and service are delivered. One reason 
for this is that many goods, e.g., food and other supplies, 
and most services, e.g., getting a haircut or going to a 
movie, are not tangible and disappear after consump- 
tion. Because they are not conserved, and also because 
they are measured in different physical units, it is not 
very practical to keep track of them. In contrast, money 
is measured in the same unit (within a given country 
with a single currency) and is conserved in local transac- 
tions ([HI) , so it is straightforward to keep track of money 
flow. It is also important to realize that an increase in 
material production does not produce an automatic in- 
crease in money supply. The agents can grow apples on 
trees, but cannot grow money on trees. Only a central 
ban k has the monopoly of chan ging the monetary base 
Mb ( McConnell and Bruel . ll996( ). (Debt and credit issues 
are discussed separately in Sec. WU ) 

Enforcement of the local conservation law ([5]) is the key 
feature for successful functioning of money. If the agents 
were permitted to "manufacture" money, they would be 
printing money and buying all goods for nothing, which 
would be a disaster. The purpose of the conservation 
law is to ensure that an agent can buy goods and prod- 
ucts from the society (the other agents) only if he or she 
contributes something useful to the society and receives 
money payment for these contributions. Thus, money 
is an accounting device, and, indeed, all accounting sys- 
tems are based on the conservation law ([5|) . The physical 
medium of money is not essential as long as the local con- 
servation law is enforced. The days of gold standard are 
long gone, so money today is truly the fiat money, de- 
clared to be money by the central bank. Money may be 
in the form of paper currency, but today it is more of- 
ten represented by digits on computerized accounts. So, 



money is nothing but bits of information, and monetary 
system represents an informational layer of the economy. 
Monetary system interacts with physical system (produc- 
tion and consumption of material goods), but the two 
layers cannot be transformed into each other because of 
their different nature. 

Unlike, ordinary economic agents, a central bank or a 
central government can inject money into the economy, 
thus changing the total amount of money in the system. 
This process is analogous to an influx of energy into a sys- 
tem from external sources, e.g., the Earth receives energy 
from the Sun. Dealing with these situations, physicists 
start with an idealization of a closed system in thermal 
equilibrium and then generalize to an open system sub- 
ject to an energy flux. As long as the rate of money 
influx from central sources is slow compared with relax- 
ation processes in the economy and does not cause hy- 
perinflation, the system is in quasi-stationary statistical 
equilibrium with slowly changing parameters. This sit- 
uation is analogous to heating a kettle on a gas stove 
slowly, where the kettle has a well-defined, but slowly in- 
creasing, temperature at any moment of time. A flux of 
money may be also produced by international transfers 
across the boundaries of a country. This process involves 
complicated issues of multiple currenci es in the world and 
their exchange rates (|McCaulevl . l2008l) . Here we consider 
an idealization of a closed economy for a single country 
with a single currency. 

Another potential problem with conservation of money 
is debt. This issue will be discussed in Sec. IVII As a 
starting point. lOragulescu and Yakovenkol (|200Cll ) consid- 
ered simple models, where debt is not permitted, which 
is also a common idealization in some econom ic liter- 
ature ()Kivotaki and Wrightl . Il993l: iMohcol . |2006[ ). This 
means that money balances of the agents cannot go be- 
low zero: m^ > for all i. Transaction ([5]) takes place 
only when an agent has enough money to pay the price: 
rrii > Am, otherwise the transaction does not take place. 
If an agent spends all money, the balance drops to zero 
rrii — 0, so the agent cannot buy any goods from other 
agents. However, this agent can still receive money from 
other agents for delivering goods or services to them. In 
real life, money balance dropping to zero is not at all 
unusual for people who live from paycheck to paycheck. 

Macroeconomic monetary policy issues, such as mone y 
supply and money demand ( Friedman and Hahnl . 119901) , 
are outside of the scope of this Chapter. Our goal is to 
investigate the probability distribution of money among 
economic agents. For this purpose, it is appropriate 
to make the simplifying macroeconomic idealizations, as 
described above, in order to ensure overall stability of 
the system and existence of statistical equilibrium. The 
concept of "equilibrium" is a very common idealization 
in economic literature, even though the real economies 
might never be in equilibrium. Here we extend this con- 
cept to a statistical equilibrium, which is characterized by 
a stationary probability distribution of money P{m), in 
contrast to a mechanical equilibrium, where the "forces" 
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FIG. 1 Histogram and points: Stationary probability distri- 
bution of money P{rn) obtained in additive random-exchange 
models. Solid curves: Fits to the exponential Boltzmann- 
Gibb s law (O . Vertical line: The initial distribution of money. 
From iDraeulescu and Yakovenkol (J2000l ) . 



of demand and supply balance each other. 



IV. THE BOLTZMANN-GIBBS DISTRIBUTION OF 
MONEY 

Having recognized the principle of local money con- 
servation, iDragulescu and Yakovenkol ()2000l ) argued that 
the stationary distribution of money P{m) should be 
given by the exponential Boltzmann-Gibbs function anal- 
ogous to Eq. (dl 



P{m) 



-m/Trr, 



(7) 



Here c is a normalizing constant, and T„i is the "money 
temperature" , which is equal to the average amount of 
money per agent: T = (m) = M/N, where M is the total 

money, and N is the number of agents.^ 

To verify this conjecture, IDragulescu and Yakovenkol 
(|200Cll ) performed agent-based computer simulations of 
money transfers between agents. Initially all agents were 
given the same amount of money, say, SlOOO. Then, a 
pair of agents (i, j) was randomly selected, the amount 
Am was transferred from one agent to another, and the 
process was repeated many times. Time evolution of 
the probability distribution of mon ey P(m) is illustrated 
i n com puter a nimation video s by IChen and Yakovenkol 
(|2007l ) and by IWrightl (|2007l) . After a transitory pe- 
riod, money distribution converges to the stationary form 
shown in Fig. [TJ As expected, the distribution is well fit- 
ted by the exponential function ([7|). 



^ Because debt is not permitted i n this model, we have M = M^, 
where Mi, is the monetary base l lMcConnell and Bruel . 0-9961) . 



Several different rules f o r Am were considered by 
iDragulescu and Yakovenkol ()2000[ ). In one model, the 
transferred amount was fixed to a constant Am = $1. 
Economically, it means that all agents were selling their 
product s for the same price Am — $1. Computer ani- 
mation ( Chen and Yakovenkol 120071 ) shows that the ini- 
tial distribution of money first broadens to a symmet- 
ric Gaussian curve, characteristic for a diffusion pro- 
cess. Then, the distribution starts to pile up around 
the 771 = state, which acts as the impenetrable bound- 
ary, because of the imposed condition ?ti > 0. As a re- 
sult, P{rn) becomes skewed (asymmetric) and eventu- 
ally reaches the stationary exponential shape, as shown 
in Fig. [H The boundary at m = is analogous to 
the ground-state energy in statistical physics. Without 
this boundary condition, the probability distribution of 
money woul d not reach a stationary st a te. Coniputer 
animations (' Chen and Yakovenkol 120071 : IWrightl . l2007l ) 
also show how the entropy of money distribution, defined 
as S/N = —^^P(mfc) In P(mfe), grows from the initial 
value S — 0, where all agents have the same money, to 
the maximal value at the statistical equilibrium. 

While the model with Am = 1 is very simple and 
instructive, it is not realistic, because all prices are 
taken to be the same. In a nothe r model considered by 
IDragules cu and Yakovenkol (2000), Am in each trans- 
action is taken to be a random fraction of the aver- 
age amount of money per agent, i.e.. Am = i'{M/N), 
where i^ is a uniformly distributed random number be- 
tween and 1. The random distribution of Am is sup- 
posed to represent the wide variety of prices for differ- 
ent products in the real economy. Computer simulation 
of t his model produces the sam e stationary distribution 
dll). IDragulescu and Yakovenkol (|200 0) also considered a 
model with firms, which hire agents to produce and sell 
products. This process results in a many-body trans- 
fer of money, as opposed to pairwise transfer discussed 
above. Computer simulation of this model also generates 
the same exponential distribution ([7|). 

Th ese ideas were further developed by Garibaldi et al\ 
(|2007) : Scalas et at (2006 ). The Bohzmann distribution 
was indepe ndently applied to social sciences by the ph ysi- 
cist JiirgenlMimtei (|2000[ ): lMimkes and Willii (|2005[) us- 
ing the Lagrange principle of maximization with con- 
straints. The exponential distrib ution of money was also 
found by the economist Martin IShubikI (|l999[ l using a 
Markov chain appro ach to strategic ma rket games. A 
long time ago, Benoit iMandelbroti (|l960l p 83) observed: 



"There is a great temptation to consider the 
exchanges of money which occur in economic 
interaction as analogous to the exchanges of 
energy which occur in physical shocks be- 
tween gas molecules." 

He realized that this process should result in the expo- 
nential distribution, by analogy with the barometric dis- 
tribution of density in the atmosphere. However, he dis- 
carded this idea, because it does not produce the Pareto 
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FIG. 2 Histogram: Stationary probability distribution of 
money in the multiplicative random-exchange model ((5)) for 
7 — 1/3. Solid c urve: The exponential Boltzmann-Gibbs law. 
From iDragulescu and Yakovcnka (|2000i ) . 



power law ( Paretol Il897l ). and proceeded to study the 
stable Levy distributions. Ironicall y, the actual economic 
data ( Yakovenko and Rosseij . |2009() do show the expo- 
nential distribution for the majority of the population. 
Moreover, the data have a finite variance, so the stable 
Levy distributions are not applicable because of their in- 
finite variance. 



V. PROPORTIONAL MONEY TRANSFERS AND 
SAVING PROPENSITY 

In the models of money transfer discussed in Sec. 
ilVi the transferred amount Am is typically indepen- 
dent of the money balances of the agents involved. A 
different model was introduced earlier b y the physicists 
llspolatov. Krapivskv. and Redneij (|l998l ) and called the 
multiplicative asset exchange model. This model also sat- 
isfies the conservation law, but the transferred amount of 
money is a fixed fraction 7 of the payer's money in Eq. 



Am, — 'ynii. 



(8) 



The stationary distribution of money in this model, com- 
pared in Fig. 12] with an exponential function, is similar, 
albeit not exactly equal, to the Gamma distribution: 



P{m)=cm^e~"'/^. 



(9) 



Eq. ^ differs from Eq. ([7]) by the power-law prefactor 
m^ . For /3 > 0, the population with low money balances 
is reduced, and i-'(O) = 0, as shown in Fig. [5] 

Essentially the same model ( Luxi . l2005i) . called the in- 
equality process, has been intr oduced and s t udied much 
earlier by the sociologist John lAngld (|l986l Il992l Il993l 



Il996l I2OO2I l2006l ). lAngld (|l986l ) associated the propor- 
tionality rule dHt with the surplus theory of social strat- 
ification (JEngelsl . Il972[ ) , which argues that inequality in 
human society develops when people can produce more 
than necessary for minimal subsistence. This additional 
wealth (surplus) can be transferred from original pro- 
ducers to other people, thus generating inequality. Angle 
found a Gamma-like distribution in numerical simula- 
tions of his m odels. Independently, the economist Miguel 
iMolicol (|2006[ ) studied conservative exchange models ^ 
where agents bargain over prices in their transactions. 
He also found a stationary Gamma-like distribution of 
money in numerical simulations. 

Another m odel with an element of propor ti onalit y was 
proposed by IChakraborti and Chakrabartil (2000). In 



this model, the agents set aside (save) some fraction of 
their money Xrrii, whereas the rest of their money bal- 
ance (1 — X)mi becomes available for random exchanges. 
Thus, the rule of exchange ^ becomes 

m'^ — Xrrii + ^(1 — X){mi + rrij), 

m'j = Xnij + (1 - 0(1 - A)(m, + ruj). (10) 

Here the coefficient A is called the saving propen- 
sity, and the random variable f is uniformly dis- 
tributed between and 1. Compu ter simulations by 
IChakraborti and Chakrabartil (|2000l ) of the model ^ 
found a stationary distribution close to the Gamma dis- 
tribution (|9|). With A 7^ 0, agents always keep some 
money, so their balances never drop to zero, thus -P(O) = 
0. 

In the subsequent papers by the Kolkata school and 
related papers, the case of random saving propensity was 
studied. In these models, the agents are assigned random 
parameter s A drawn from a uniform distribution b etwee n 
and 1 (|Chatteriee. Chakrabarti. and Manual . |2004|) . 
It was found that this model produces a power-law tail 
P(m) oc Ijim? at high m. The reasons for stability of 
this law wer e understood using the Boltzmann kinetic 
equat ion (IChatteriee. Chakrabarti. and Stinchcombe 



2005: 



Das a nd Yarlagaddal 



Repetowicz. Hutzler. an d Richmondl, 

but most elegantly in the mean-fiel d 

(iBhattacharvva. Chatteriee^ and Chakrabartil 

IChatterjee and Chakrabarti 120071 : iMohantvl 

The fat tail originates from the agents whose saving 

propensity is clo se to 1, who hoard money and do not 

give i t back (IPatriarca. Chakraborti. and Germanoj . 




I2OO6I : IPatriarca e^aLl . 120051 )" 
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ous mathema tical treatment of . 

given by Diiring. Matthes. and Toscanil (120081): 



more rigor- 
prq blem was 



(|2007l ): iMatthes and Toscanil 



iDiiring and Toscani 
(l2008f ). 

As a further extension, iDragulescu and Yakovenkoj 
(2000) considered a model with taxation, which also has 
an element of proportionality. The Gamma distribution 
was also studied for con servative mode l s with in a simple 
Boltzmann approach bv iFerrerol ( 2004120051) and, using 



more complicate d rules of exchange motivate d by po lit- 
ical economy, by IScafetta. Picozzi. and West! (|2004al |bl) . 
Another extension of these studies includes not only 
money transfers, but also tra nsfers of a commod- 
ity, for which money is paid (JAusloos and Pekalski , 
20071 Chakraborti. Pradhan. and Chakrabarti . 



200li I Chatteriee and Chakrabartil 120061: iLuxl. 12009: 
Silver. S lud. and Takamoto, 200 ^. For a more detailed 
review of these models, see lYakovenko and Rosseij 
(l2009h 



The stationary distribution of money ([9]) is different 
from the simple exponential formula ([7|). The origin of 
this difference c an be understood from the Boltzmann 
kinet ic equation (JLifshitz and Pitaevskii Il98ll : IWannien . 
Il987l ). This equation describes time evolution of the dis- 
tribution function P{m) due to pairwise interactions: 



dP{m) 

dr~ 



{-/[: 



rti.ni \—>\7n — l\.7n 



-A]P(m)P(m') (11) 



+ /[m--A,m'+A]^[m,m']^(™ " ^)P{m' + A)} dm' dA. 

Here /[m,m']-i-[m-A,m'+A] IS the probability of transfer- 
ring money A from an agent with money m to an agent 
with money m' per unit time. This probability, multi- 
plied by the occupation numbers P{m) and P[m')^ gives 
the rate of transitions from the state [m, m'] to the state 
[m — A, m' -I- A]. The first term in Eq. (|TT|) gives the de- 
population rate of the state m. The second term in Eq. 
(fTTjl describes the reversed process, where the occupation 
number P(m) increases. When the two terms are equal, 
the direct and reversed transitions cancel each other sta- 
tistically, and the probability distribution is stationary: 
dP{m)/dt = 0. This is the principle of detailed balance. 
In physics, the fundamental microscopic equations of 
motion obey the time-reversal symmetry. This means 
that the probabilities of the direct and reversed processes 
are exactly equal: 



/[n 



>[m — A,; 



-A] — /[m-A,m'+A]->[m,m']- (12) 



When Eq. (fT^ is satisfied, the detailed balance condi- 
tion for Eq. (fTTjl reduces to the equation P{m)P{m') — 
P{m — A)P{m' + A), because the factors / cancels out. 
The only solution of this equation is the exponential func- 
tion P(rn) = cexp(— 77i/Tm), so the Boltzmann-Gibbs 
distribution is the stationary solution of the Boltzmann 
kinetic equation (fTTjl . Notice that the transition prob- 
abilities (fT2j) are determined by the dynamical rules of 
the model, but the equilibrium Boltzmann-Gibbs distri- 
bution does not depend on the dynamical rules at all. 
This is the origin of the universality of the Boltzmann- 
Gibbs distribution. We see that it is possible to find the 
stationary distribution without knowing details of the dy- 
namical rules (which are rarely known very well) , as long 
as the symmetry condition (|12p is satisfied. 

The models considered in Sec. IIVI have the time- 
reversal symmetry. The model with the fixed money 
transfer A has equal probabilities (fT^ of transferring 
money from an agent with the balance m to an agent 



with the balance m' and vice versa. This is also true 
when A is random, as long as the probability distribution 
of A is independent of m and m'. Thus, the stationary 
distribution P[m) is always exponential in these models. 
On the other hand, in the model ([8]), the time-reversal 
symmetry is broken. Indeed, when an agent i gives a 
fixed fraction 7 of his money mi to an agent with bal- 
ance rrij, their balances become (1 — 7)77^^ and mj +"fmi. 
If we try to reverse this process and appoint the agent j 
to be the payer and to give the fraction 7 of her money, 
jirUj + ^rrii), to the agent i, the system does not return 
to the origin al configuration [mi^rrij]. As emphasized by 
lAngld (|2006() . the payer pays a deterministic fraction of 
his money, but the receiver receives a random amount 
from a random agent, so their roles are not interchange- 
able. Because the proportional rule typically violates 
the time-reversal symmetry, the stationary distribution 
P{m) in multiplicative models is not exponential. 

These examples show that the Boltzmann-Gibbs dis- 
tribution does not necessarily hold for any conserva- 
tive model. However, it is universal in a limited sense 
for a broad class of models that have time-reversal 
symmetry. In the absence of detailed knowledge of 
real microscopic dynamics of economic exchanges, the 
semiuniversal Boltzmann-Gibbs distribution ([7]) is a 
natural starting point. Moreove r, the assumption of 
iDragulescu and Yakovenkq ([2000) that agents pay the 
same prices Am for the same products, independent of 
their money balances m, seems very appropriate for the 
modern anonymous economy, especially for purchases 
over the Internet. There is no particular empirical ev- 
idence for the proportional rules © or ([TU]). By further 
modifying the rules of money transfer and introducing 
more parameters in the models, it i s possible to obtain 
even more complicated distributions (jScafetta and Westi . 
I2OO7I ). However, parsimony is the virtue of a good math- 
ematical model, not the abundance of additional assump- 
tions and parameters, whose correspondence to reality is 
hard to verify. 



VI. MODELS WITH DEBT 

Now let us discuss how the results change when debt is 
permitted.^ From the standpoint of individual economic 
agents, debt may be considered as negative money. When 
an agent borrows money from a bank (considered here as 
a big reservoir of money), '^ the cash balance of the agent 



^ The i deas presented here are quite similar to those by ISoddvl 
lll933h . Frederick Soddy, the Nobel Prize winner in chemistry 
for his work on radioactivity, argued that the real wealth is de- 
rived from the energy use in transforming raw materials into 
goods and services, and not from monetary transactions. He 
also warned about dangers of excessive debt and related "virtual 
wealth" resulting in the Great Depression. 

^ Here we treat the bank as being outside of the system consisting 
of ordinary agents, because we are interested in money distri- 



(positive money) increases, but the agent also acquires 
a debt obligation (negative money), so the total balance 
(net worth) of the agent remains the same. Thus, the act 
of borrowing money still satisfies a generalized conserva- 
tion law of the total money (net worth), which is now 
defined as the algebraic sum of positive (cash Af ) and 
negative (debt D) contributions: M — D = Mb, where Mb 
is the origi nal amount of money in the system, the mon- 
etary base ([McConnell and Brud . Il996( ). After spending 
some cash in pairwise transactions ([5]) , the agent still has 
the debt obligation (negative money) , so the total money 
balance to^ of the agent (net worth) becomes negative. 
We see that the boundary condition m^ > 0, discussed in 
Sec. lIIH does not apply when debt is permitted, so m = 
is not the ground state any more. The consequence of 
permitting debt is not a violation of the conservation law 
(which is still preserved in the generalized form for net 
worth) , but a modification of the boundary condition by 
permitting agents to have negative balances ?Tii < of net 
worth. A more detailed discussion of positive and neg- 
ative money and the book-keeping accounting from the 
econo physics point o f vie w was presented by the phys icist 
Dieter lBraunr(|200l[ ) and lFischer and BraunI (|2003allbl ). 

Now we can repeat the simulation described in Sec. 
IIVI without the boundary condition to > by allowing 
agents to go into debt. When an agent needs to buy a 
product at a price Am exceeding his money balance rrii, 
the agent is now permitted to borrow the difference from 
a bank and, thus, to buy the product. As a result of 
this transaction, the new balance of the agent becomes 
negative: m.^ = rrii — Am < 0. Notice that the local con- 
servation law ([5]) and 1^ is still satisfied, but it involves 
negative values of m. If the simulation is continued fur- 
ther without any restrictions on the debt of the agents, 
the probability distribution of money P(m) never stabi- 
lizes, and the system never reaches a stationary state. As 
time goes on, P{m,) keeps spreading in a Gaussian man- 
ner unlimitedly toward to. = -|-cx) and to — —00. Because 
of the generalized conservation law discussed above, the 
first moment (m) = Mb/N of the algebraically defined 
money m, remains constant. It means that some agents 
become richer with positive balances m > at the ex- 
pense of other agents going further into debt with nega- 
tive balances to < 0, so that M = Mb + D. 

Common sense, as well as the experience with the cur- 
rent financial crisis, tells us that an economic system can- 
not be stable if unlimited debt is permitted.* In this case, 
agents can buy any goods without producing anything in 
exchange by simply going into unlimited debt. Arguably, 
the current financial crisis was caused by the enormous 



bution among these agents. The debt of agents is an asset for 
the bank, and deposits of cash into the bank are liabilities of 
the bank JMcConnell and Bruel 119961 ) . We do not go into these 

details in order to keep our presentation simple. 

* In qu alitative agreement with the conclusions by iMcCaulevI 
l l2008f) . 



N=500, M=5*10^, time=4*10^. 




4000 
Money, m 

FIG. 3 Histograms: Stationary distributions of money with 
and without debt. The debt is limited to md = 800. 
Solid curves: Fits to the Boltzmann-Gibbs laws with the 
"money temperatures" Tm = 18 00 and r,„ — 1000. From 
iDraeulescu and Yakovenkg (J2000l '). 
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FIG. 4 The stationary distribution of money for the required 
reserve ratio R — 0.8. The distribution is exponential for both 
positive and negative money with different "temperatures" T+ 
and T-, as illustrated by t he inset on log-linear scale. From 
IXi. Ding, and Wand l|2005l ). 



debt accumulation in the system, triggered by subprime 
mortgages and financial derivatives based on them. A 
widely expressed opinion is that the current crisis is not 
the problem of liquidity, i.e., a temporary difficulty in 
cash flow, but the problem of insolvency, i.e., the inher- 
ent inability of many participants to pay back their debts. 

Detailed discussion of the current economic situation 
is not a subject of this paper. Going back to the idealized 
model of money transfers, one would need to impose some 
sort of modified boundary conditions in order to prevent 
unlimited growth of debt a nd to ensure o verall stability 
of the system. iDragulescu and Yakovenko (2000 ) con- 
sidered a simple model where the maximal debt of each 
agent is limited to a certain amount rud- This means 
that the boundary condition tti^ > is now replaced by 
the condition m^ > —rud for all agents i. Setting inter- 
est rates on borrowed money to be zero for simplicity, 
IDragulescu and Yakovenkg (J200Cl() performed computer 
simulations of the models described in Sec. IIVI with the 
new boundary condition. The results are shown in Fig. 
[3l Not surprisingly, the stationary money distribution 
again has the exponential shape, but now with the new 
boundary condition at m = —md and the higher money 
temperature Td = rud + Mb/N. By allowing agents to go 
into debt up to nid, we effectively increase the amount of 
money available to each agent by nid- So, the money tem- 
perature, which is equal to the average amount of effec- 
tively available money per agent, increases correspond- 
ingly^ 

IXi. Ding, and Wang] ( 20051 ) considered another, more 
realistic boundary condition, where a constraint is im- 
posed not on the individual debt of each agent, but 
on the total debt of all agents in the system. This is 
accomplished via the required reserve ratio _R, which 



is briefly explained below ( McConnell and Bruel Il996() . 
Banks are required by law to set aside a fraction R of 
the money deposited into bank accounts, whereas the 
remaining fraction 1 — R can be loaned further. If the 
initial amount of money in the system (the money base) 
is Mh, then, with repeated loans and borrowing, the to- 
tal amount of positive money available to the agents in- 
creases to M = Mh/ R, where the factor 1/R is called 
the money multiplier (JMcConnell and Bruel [l996( ) . This 
is how "banks create money". Where does this extra 
money come from? It comes from the increase in the to- 
tal debt in the system. The maximal total debt is given 
by 13 = Mb/R — Mb and is limited by the factor R. When 
the debt is maximal, the total amounts of positive, Mb/ R, 
and negative, Mb{l — R)/R, money circulate among the 
agents in the system, so there are two constraint s in the 
model considered bv lXi. Ding, and Wand (|2005( ). Thus, 
we expect to see the exponential distributions of positive 
and negative money characterized by two different tem- 
peratures: r+ = Mb/RN and T_ = Mb{l - R)/RN. 
This is exactly what was found in computer simula- 
tions by IXi. Ding, and Wangj ( 20051 ). as shown in Fig. 
H Similar two-sided dist ributions were also found by 



Fischer and BraunI (|2003aD 



However, in the real economy, the reserve requirement 
is not effective in stabilizing total debt in the system, be- 
cause it applies only to deposits from general public, but 
not from corporations. Moreover, there are alternative 
instruments of debt, including derivatives and various un- 
regulated "financial innovations" . As a result, the total 
debt is not limited in practice and can potentially reach 
catastrophic proportions. Here we briefiy discuss several 
models with non-stationary debt. Thus far, we did not 
consider the interest rates. IDragulescu and Yakovenkol 
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(2000) studied a simple model with different interest rates 
for deposits into and loans from a bank. Computer simu- 
lations found that money distribution among the agents 
is still exponential, but the money temperature slowly 
changes in time. Depending on the choice of param- 
eters, the total amount of money in circulation either 
increases or decreases in time. Interest amplifies desta- 
bilizing effect of debt, because positive balances become 
even more positive and negative even more negative due 
to accruement of interest. A more sophisticated m acroe- 
conoin i c mo del studied by the economist Steve iKeenI 
( 19951 |2000( ) exhibits debt- induced breakdown, where 
all economic activity stops under the burden of heavy 
debt and cannot be restarted without a "debt morato- 
rium" . The interest rates were fi xed in these models and 
not a djusted self-consistently. ICockshott and Cottrelll 
( 20081 ) proposed a mechanism, where the interest rates 
are set to cover probabilistic withdrawals of deposits 
from a bank. In an agent-base d simulation of the model, 
ICockshott and Cottrelll (|2008l ) found that money supply 
first increases up to a certain limit, and then the econ- 
omy experiences a spectacular crash under the weight of 
accumulated debt. 

In the absence of a nominal limit on maximal debt, 
bankruptcy provides a mechanism for debt stabilization. 
When the debt of an agent becomes too large, the agent 
will not be able to borrow any more money and will not 
be able to pay the debt back, so he or she will have 
to declare bankruptcy. Bankruptcy erases the debt of 
the agent (the negative money) and resets the balance 
to zero. However, somebody else (a bank or a lender) 
counted this debt as a positive asset, which also be- 
comes erased. In the language of physics, creation of 
debt is analogous to particle-antiparticle generation (cre- 
ation of positive and negative money), whereas cancella- 
tion of debt by repayment or by bankruptcy corresponds 
to particle-antiparticle annihilation (annihilation of posi- 
tive and negative money) . The former process dominates 
during economic bubbles (booms) and represents mone- 
tary expansion, whereas the latter dominates during the 
subsequent recessions (busts) and represents monetary 
contraction. Bankruptcy is the crucial mechanism for 
stabilizing money distribution, but it is often overlooked 
by the economists. Interest rates are meaningless without 
a mechanism specifying when bankruptcy is triggered. 

After lending money out, the lender has the burden 
of collecting the debt from the debtor. Thus, the act of 
lending creates a persistent link (a string) between the 
lender and the borrower, as emphasized in the Chapter 
by Heiner Ganssmann in this Volume. This is in contrast 
to payments by positive money for goods and services, 
which are final and do not leave any persistent link be- 
tween the agents after the transaction. Invention of the 
infamous collateralized debt obligations (CDO) obscured 
connections between lenders and borrowers by randomiz- 
ing and anonymizing their pools. It destabilized the sys- 
tem by inviting unsustainable debt and made bankruptcy 
proceedings extremely difficult because of the scrambled 



identities of lenders and borrowers. 

As argued above, boundary conditions are crucial for 
stabilizing money distribution. Typically, a lower bound 
is imposed, but not an upper bound (in a capitalist, as 
opposed to a socialist, society)."" This asymmetry is very 
important for stability of a monetary system. Numer- 
ous attempts were made to create alternative community 
money from scratch and most of them failed. In such a 
system, when an agent provides goods or services to an- 
other agent, their accounts are credited with positive and 
negative tokens, as in Eq. ([S])- However, because the ini- 
tial global money balance is zero in this case, the prob- 
ability distribution of money P{m) is symmetric with 
respect to positive and negative m. Unless a boundary 
condition is imposed on the lower side, P{m) will never 
stabilize. Some agents will accumulate unlimited nega- 
tive balance by consuming goods and services and not 
contributing anything in return, thus undermining the 
system. In contrast, when a central government creates 
positive money by fiat and forces its usage by demanding 
that taxes are paid with this money, it creates a viable 
monetary system, as discussed in the Chapter by Ran- 
dall Wray in this Volume. Thus, taxation is an essential 
ingredient for vitality of a monetary system. 



VII. CONCLUSIONS AND PERSPECTIVES 

In this Chapter, we have demonstrated that random 
transfers of money in economic transactions between oth- 
erwise equal economic agents produce a broad and highly 
unequal probability distribution of money among the 
agents. In additive models, the probability distribution 
of money is exponential and similar to the Boltzmann- 
Gibbs distribution of energy in statistical physics. Multi- 
plicative models produce a Gamma-like distribution and 
a power-law tail for random saving propensity. Local 
conservation of money in transactions between agents is 
crucial for the accounting function of money. Ordinary 
economic agents can only receive and give money, but 
cannot produce it. Boundary conditions are necessary 
in order to achieve a stable probability distribution of 
money. Without debt, zero money balance is the bound- 
ary. When debt is permitted, some sort of restriction 
on the negative money balances must be imposed, either 
by limiting individual or collective debt or by setting up 
conditions for bankruptcy. When debt is unlimited, the 
system is unstable and does not have a stationary state. 

It would be very interesting to compare these the- 
oretical conclusions with empirical data on money 
distribution. Unfortunately, it is very difficult to obtain 



^ If an upper limit is imposed instead of a lower limit, the money 
temperature in Eq. ([7]l becomes negative Tm < 0, so the slope in 
Fig. [T] changes to dP/dm > 0, which is known in physics as the 
inverse popula tion. The case with both upper a nd lower limits 



was studied bv lDragulescu and Yakovenkd ll200Cll ) 



r, 



such data. The probability distribution of balances on 
deposit accounts in a big enough bank would be a rea- 
sonable approximation for money distribution among the 
population. However, such data are not publicly avail- 
able. In contrast, plenty of data are available on income 
distribution from the tax a gencies. Quantitative analysis 
of such data for the USA ( Yakovenko and Rosseii |2009[ ) 
shows that the population consists of two distinct social 
classes. Income distribution follows the exponential law 
for the lower class (about 97% of population) and the 
power law for the upper class (about 3% of population) . 
Although social classes have been known since Karl 
Marx, it is interesting that they can be recognized by 
fitting the empirical data with simple mathematical 
functions. Using sophisticated models of interacting 
econo r nic ag ents, the computer scientist Ian IWrighi 
( 20051 I2OO9I) demonstrated emergence of two classes 
in agent-based simulations of initially equal agents. 
This work has been further developed in the book by 
ICockshott. Cottrell. Michaelson. Wright, and Yakovenkol 



()2009l ). integrating economics, computer science, and 
physics. 

Nowadays, money is typically represented by data bits 
on computer accounts, which constitute the informa- 
tional layer of the economy. In contrast, material stan- 
dards of living are determined by the physical layer of 
the economy. In the modern society, physical standards 
of living are largely determined by the level of energy 
consumption and are widely di fferent around the globe. 
iBaneriee and Yakovenkol ( 2010( ) found that the probabil- 
ity distribution of energy consumption per capita around 
the world approximately follows the exponential law. So, 
it is likely that the energy consumption inequality is gov- 
erned by the same principles as the money inequality. 
The energy/ecology and financial/economic crises are the 
biggest challenges faced by the mankind today. There is 
an urgent need to find ways for a manageable and realis- 
tic transition from the current breakneck growth-oriented 
economy, powered by ever-expanding use of fossil energy 
fuels, to a stable and sustainable society, based on renew- 
able energy and balance with the Nature. Undoubtedly, 
both money and energy will be the key factors shaping 
up the future of human civilization. 
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